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Abstract. - The homogeneous entropy for continuous systems in nonextensive statistics reads 
S" =k B (l- (K faTp 1/q (T)) q )/{l -q), where T is the phase space variable. Optimization of S„ 
combined with normalization and energy constraints gives an implicit expression of the distribution 
function p(r) which can be computed explicitly when the Hamiltonian reduces to its kinetic part. 
We examine the meaning of the q-ideal gas and we compute its properties such as the energy 
fluctuations and the specific heat. Similar results are also presented using the formulation based 
on the Tsallis entropy. From the analysis, we find that the validity of the nonextensive formalism 
for the q-ideal gas is restricted to the range q < 1, which raises the question of the formal validity 
range for continuous Hamiltonian systems. 



Introduction. — Many phenomena in natural systems and in laboratory experiments 
are observed and measured under non-equilibrium conditions, and therefore do not obey the 
standard statistical mechanics description. In particular the distributions which characterize 
such systems are not Boltzmann-likc and do not follow from linear response. Instead these 
distributions exhibit "fat tails" and power law decays and often they can be fitted by q- 
cxponcntial functions which generalize the usual Boltzmann exponential distribution [1] . A 
recurring question is what deviation from standard statistical mechanics gives rise to this 
behavior, which amounts to the question of the emergence of "statistics from dynamics" as 
emphasized by E.G.D. Cohen [2]. There are several possible analytical developments from 
which q-cxponcntial distributions can be obtained: superstatsitics [5] by statistical average 
over the %-square distribution of an intensive variable, nonlinear response theory [3] by 
the solution of the generalized Fokker-Planck equation, and nonextensive statistics [5] by 
optimization of the generalized entropy. 

It was precisely the original idea of nonextensive statistics introduced by Tsallis about 
20 years ago [5] to develop a statistical mechanical theory for systems out of equilibrium 
where the Boltzmann distribution no longer holds, and to generalize the Boltzmann entropy 
by a more general function S q while maintaining the formalism of thermodynamics. From 
a practical viewpoint, the nonextensive statistics formulation appeared to be of interest 
because maximization of the generalized entropy under the usual constraints (normalized 
probabilities, fixed internal energy) yields the q-cxponcntial distribution which has been 
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successfully used to describe distributions observed in a large class of phenomena [6] . Indeed 
for a certain range of values of the index q, these g-exponential distributions exhibit a power 
law decay (when q > 1), a feature observed in a large class of experimental phenomena 
which cannot be straightforwardly interpreted in the context of classical theories. At the 
same time, a large literature concerning the internal self-consistency of the nonextensive 
formalism has developed addressing such questions as the stability of the entropy functional 
[?J, the method of calculating averages [10] and the positivity of the specific heat [8lfTT] . 

Many explicit applications of the nonextensive formalism involve the assumption of in- 
dependent particles: e.g., non- interacting particles that can occupy a set of discrete energy 
levels or, in the classical case, the ideal gas. These applications may appear paradoxical as 
the assumption of non-extensivity implies an interaction between components of the system. 
An alternative point of view, adopted here, is that there is no paradox because the underlying 
physical system does involve interactions (even long-ranged interactions). The independent 
particles are not the constituents of the physical system, but rather are understood to be 
quasi-particles in which the effect of the interactions are, to a first approximation, included 
in their properties (effective mass, statistics, etc.). With this point of view, the adoption of 
the nonextensive formalism is another part of the effective one-body description required to 
account for aspects of the interactions that cannot be otherwise modeled. Accordingly we 
revisit the formalism for the q- ideal gas. 

Even this simple case of the g-ideal gas has been subject to questions of internal consis- 
tency. In particular, that the nonextensive formalism has limited range of validity [10] . gives 
negative specific heats [TT] and even, recently, negative values for the second cumulant of 
the energy (a positively defined quantity!) [8]. Here we re-examine these issues using both 
the Tsallis entropy and, as suggested recently [7J, the homogeneous entropy. 

The homogeneous entropy. — We start with the homogeneous entropy (or normal- 
ized Tsallis entropy) Sjf which was proven to be stable against small perturbations in the 
probability distribution function p(T) while the Tsallis entropy S q is not [7]. For continuous 
systems the f/-entropy is given by 



s?=k B - ^ yLn , a) 



i-q 

where q is the index characterizing the entropy functional, T denotes the phase space vari- 
able, and K must be a quantity with the dimensions of [r] q , i.e. K = h K a ) with N, 
the number of degrees of freedom of the system with dimension D and Hamiltonian H. In 
the limit q — > 1, the classical Boltzmann-Gibbs formulation is retrieved^ 

Optimization of the ii-entropy |T]) with the normalization and energy constraints: 

1 = J p{T)dT ■ U = J p (r) HdT , (2) 

by the method of Lagrange multipliers leads to 



JW l 



?f -a(j p(r)dT-l)-pn p(T)HdT-U 



-q(K$p^{T)dT) q - l Ky^(T) 



which is solved to give 



-a-fiH, (3) 



p(r) 



(q-l)a {q-l)P 




\{K j p s -/i{Y)dT) q X K (K j p 1 / '9(T) dT) q 1 K 
Z\ {a'+P'H)? 5 , (4) 



1 For simplicity the Boltzmann factor ks will be omitted and reincluded explicitly when necessary. 
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where Z q = K J p 1 ^ (T) dT, and a' = (q- l)a/K and /3' = (q - 1)P/K. 

p (r) is the physical probability distribution, which must be real, positive and normaliz- 
able; so we must have 

1 = 2% [ (a' + (3'H)f^dT. (5) 



We will also consider the first moments of the Hamiltonian 

(H m ) = J p(T) H m dT = Z\ J {a + /3'H)f* H m dT ; (6) 
in particular we are interested in 

J {a' + (3'H) — dT J (a' + j3'H) ~ dT 

where we used the normalization condition in ([2]). So the integrals to be considered have 
the form 

I m = I (a' +pH)f* H m dT, (8) 



and the result of the integration will depend on the sign of the Lagrange multipliers; therefore 
we must consider the following possible cases: (i) a' > and j3' < 0, (ii) a' > and j3' > 0, 
and (hi) a' < and /J' > (If both a',/3' < 0, there is no solution to the normalization 
condition ([5]))- 

The q-ideal gas. — Besides the physical meaning of the q-ideal gas which was explained 
in the introductory section, it is legitimate to discuss its validity in the context of the 
nonextensive formalism because, if the formalism is to be used for continuous Hamiltonian 
systems, it should first pass the test of the g-ideal gas (as in classical statistical mechanics). 

For the g-ideal gas, the Hamiltonian reduces to its kinetic part and the configuration 
integral in ([8]) is straightforward and yields a factor given by the space volume V N . With a 
change of variable X = ^ ([SJ for a' > and (3' < (case (i)) becomes 

1 / I \ i x+ m poo 

where B(k, I) is the Beta function provided q < 1 (i.e. a < 0, and (3 > 0), and excludes the 
possibility q > 1 (with a > and /3 < 0). Then we have 



IL G ( a 



m B 



1 ND 
1-q ' 2 



(10) 



which, with ([7]), gives 



a 



and 



U-?' 2 



4 U 



2 



ND(2+(l-q) {2 + ND)) 



(12) 
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Note that (H 2 ) — (H) 2 is always positive since the q index must be q < 1. It also follows 
from these results that the explicit expression of the distribution function for the g-ideal gas 
is given by 



■■<! 
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or, with the notation exp^ = (! + (! — q)x)] r 



p (r) = f 5r ff m ,r ■ (14) 

J exp ? -g| (H - U) dT 

Noting that Z q= \ = 1, it is clear that for q = 1, one retrieves the classical exponential 
distribution. 

Proceeding along the same lines for case (ii): a',/3' > 0, i.e. a, (3 > with q > 1, or 
a, P < with q < 1, we obtain 



1 2 ' g-1 2 J 



(15) 



if and only if 1 < q < 1 + lVO | — 1 . This gives 

« = u( 2 - 1 

p U \ND(q-l) 

4 U 2 

^-^ 2 = ND(2-(q-l)(2 + ND))' (16) 

which is valid (positive definite) when 1 < q < 1 + N p +2 ■ Notice that this range of the 
q index is vanishingly small for ND >> 1 and therefore physically negligible. In this case 
a, P > 0, but the case a, P < (with q < 1) is excluded. 
For a' < and P' > (case (iii)), we have 



/I f i \ -^r-+m 
1 „^S_ / a 



7 m = F^ DiV (2m)^/ 2 ^'|^f^iJ jf (-l + JSO^Jf*- 1 ^, (17) 

which, whether g < 1 or g > 1, has no solution. So the cases g > 1 with a < 0, and P > 0, 
and g < 1 with a > 0, and /3 < are excluded. 

In summary, we have shown that, except for the physically negligible range 1 < q < 
1 + N p +2 1 the distribution function p(T) for the g-ideal gas is normalizable only for q < 1 
(with a < 0, and /? > 0), and that, contrary to some recent claim [8], the positivity of the 

energy mean squared fluctuations ({H — (H)) 2 \ is always satisfied. 



If, instead of the homogeneous entropy ([!} , we start from the Tsallis entropy [5] for con- 
tinuous systems S q = — - — - , and use the same optimization procedure ([2]) (except 

that U must then be computed with the escort average U = ^jp^l^dr ) ; we obtain the 
distribution function 

Pt (r) = (an - a/r(l - 9 ) A (H - U)j ^ , (18) 

where all — and Zt — K J p q (T) dT. Therefrom performing the computation for the 

q-idcal gas [9j leads to conclusions that are the same as above and are in essential agreement 
with some results by Abe [TUJ[TTj ; in particular we find that the normalized distribution 
function exists only for q < 1 (besides the physically vanishingly small (for N >> 1) range 
1 < q < 1 + 7715+2) w ith the additional observation that pt(F) has a singular point at q = 0. 
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Thermodynamic quantities. — We now evaluate the homogeneous entropy starting 
from (fTJ) rewritten as 

S? = k B i— 1 , (19) 

with 

Kj(a' + f3'H)t^dr) q 

Z\ = ± ^ '— = K" - . (20) 

/ (a' + p'H)\-« dT 7 ° 

For the q-ideal gas with q < 1 and f3 > 0, using ©, we find 



1-q ' 2 



where 



and 



|/3'| V (1-?)JVD 
Combining these results, we obtain 



with 



f ND _ q 



and 

3jr _ fc 1 - if 9 (V; g ) t/^" 1 ) 
1-9 

It follows that the thermodynamic temperature of the g-ideal gas is given by 



(22) 



U (1 + 7, ' )■ (23) 



9 - ^(l-,)^ ' ^ 
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(25) 



S g = fc B ; : • (26) 



1 dS% , „„„„-,„ , ND, 



k B Km H (V;q)— U— (27) 



In the limit g -> 1, A' = 1 and R H (V; <? — > 1) = 1, so that for the classical ideal gas, where 
[7 = ^P-ksT, we retrieve the expression f§ = ^ The specific heat is then readily obtained 



dT H \ 1 , ND „„ „,„ , C/^" 1 ) 



C v = (^7r) = fc B — ^j? g (F;g) i ■ - , (28) 



which is always positive for q < 1, and, for 5 = 1, gives the classical result Cy = fcs 
Note that using ([T2"]). ([23]) and ([2"7]). ([55]) can also be written as 



, fl _ ( ff 2 ) - (H ) 



2 
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with C q = Z q q f 1 + - 1 ^ ^ J , which generalizes the expression of the specific heat given 

in terms of the energy fluctuations Cy = ((AE) 2 ) / (fc^T 2 ). 

When we perform the same computation with the Tsallis formulation [9], we find the 
specific heat 

ND [/(i-s)^ 2 
C v = k B — K R(V; g) 1 ^ {1 _ q)mi , (30) 

which, in the limit q — > 1, gives the classical expression for Cy, but where the denominator 
is negative for q < 1, except when q = 1 — e with e < << 1. So, except in this narrow 
range, the Tsallis entropy formalism gives a negative specific heat for the q-ideal gas. 

Concluding comments. — We have shown that optimization of the ii-cntropy for 
continuous Hamiltonian systems combined with normalization and energy constraints gives 
an expression for the distribution function which is computed explicitly for the g-idcal gas 
and that, in the thermodynamic limit, the distribution function exists in the q < 1 index 
range. We have also shown (i) that in this range the mean squared energy fluctuations are 
always positive, in contradition to recent claims that were a result of not taking into account 
the existence of intermediate integrals in the evaluation [5] , and (ii) that in the usual Tsallis 
formulation the specific heat of the q- ideal gas is negative for q < 1. We conclude that 
the use of the non-extensive formalism to " explain" observed ^-exponential distributions on 
the basis of non-interacting quasi-particlcs is problematic when q > 1, the range where the 
g-cxponential function exhibits power law decay. Furthermore, in the range q < 1 where the 
normalized distribution function exists, the Tsallis formalism is also questionable as it gives 
a negative specific heat for the q-idcal gas. Its applicability to Hamiltonian systems with 
continuous canonical variables has also been questioned recently by Abe from a different 
viewpoint |12j . 

As discussed in the introduction, there are two aspects to the noncxtensive approach to 
the study of noncquilibrium systems, (i) Noncxtensive statistics has been applied succcsfully 
to analyze and to interpret observations in Hamiltonian systems which exhibit power law 
decay [6] ; these interpretations are based on phcnomenological analyses in accordance with 
g-exponential distributions, (ii) The nonextensive formalism was constructed on the basis 
of a few axioms and accordingly should develop with self-consistency. Our analysis suggests 
that the range of validity of the latter is limitted. It may be that there would be less 
restriction if the formalism were developed with interacting particles but in this case, the 
reason for assuming the nonextensive formalism becomes unclear. This however does not 
preclude the pragmatic application of nonextensive statistics in phenomenological analyses 
of experimental results. 
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